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Abstract 

The recent developments on the "vector manifestation" of chiral symmetry by 
Harada and Yamawaki provide a compelling evidence for, and "refine," the in-medium 
scaling of hadronic properties in dense/hot matter (call it "BR scaling") proposed by the 
authors in 1991. We reinterpret the Harada- Yamawaki result obtained in a Wilsonian 
renormalizat ion-group approach to hidden local symmetry theory matched to QCD at 
near the chiral scale in terms of the Nambu-Jona-Lasinio model and predict that the 
vector meson mass in medium should scale m*/m p ~ ((qq)* / (qq)) 1 ^ 2 from n = up to 
~ n = no (where no is nuclear matter density) and then go over to the linear scaling 
in* /m p ~ (qq)* / (qq) up to the chiral transition density n c at which the mass is to van- 
ish (in the chiral limit) according to the vector manifestation fixed point. In the regime 
of the linear scaling above n , the (vector) gauge coupling constant should fall linearly 
in (qq)*, vanishing at the critical point, with the width of the vector meson becoming 
steeply narrower. We suggest the in-medium vector meson mass, both parametric and 
pole, which goes smoothly to zero with increasing density, to be an order parameter 
for chiral symmetry restoration. Some issues related to recent RHIC observations are 
qualitatively discussed in the appendices. Our main conclusion that follows from these 
considerations is that the movement towards chiral restoration can be reliably described 
in Nambu-Jona-Lasinio mean field, with constituent quarks as variables, although there 
is probably density discontinuity following chiral restoration. 



1 Introduction 



In their "vector manifestation (VM)" of chiral symmetry, Harada and Yamawaki [1] 
have determined the parameters of their hidden local symmetry (HLS) by matching to 
QCD at a sufficiently high scale A > Arrf^ so that the HLS (hadronic) effective theory over- 
laps with QCD. When applied to hadrons in medium, masses and coupling constants in the 
hadronic theory are then "parameteric," dependent intrinsically on density/temperature. 
Renormalization group (RG) equations determine how they flow as scale is changed. In 
this note, we show that the Nambu-Jona-Lasinio (NJL) theory is a caricature theory with 
constituent quarks (or quasiquarks) as variables, which mimics most aspects of the vec- 
tor manifestation, especially the parametric behavior of the masses. NJL has the same 
quadratic divergence as the HLS, the parametric dependence of the parameters being run 
chiefly by the flow of the in-medium condensate (qq)*- Furthermore the NJL model makes 
clear the nonperturbative (solitonic) nature of the vacuum, although the fixed-point struc- 
ture of the HLS cannot be constructed until mesons are "integrated in." (For a review 
on NJL, in particular the bosonization thereof, see, e.g., Vogel and Weise [2]). The NJL 
cutoff i\njl can be interpreted as the QCD matching scale for constituent quarks. 

Introduction of vector mesons immediately raises the QCD matching scale to A ~ 
4:71 f n , the appropriate one for Wilsonian matching, in a self-consistent way. 

Our interpretation of the HLS vector manifestation of chiral symmetry clarifies: (1) 
That the use in Brown- Rho (BR) scaling [3, 4] of medium dependent masses is correct, al- 
though the approximately equal scalings such as m*/m p m f*/f n are only correct as long 
as RG flow of /* dominates, at low densities or temperatures, the flow in other variables, 
especially that of the (hidden) gauge coupling g, becoming important at higher densities 
and temperatures.; (2) Once the parameters of the Lagrangian that are intrinsically de- 
pendent on density or temperature are established, the configuration mixing such as the 
Dey-Eletsky-Ioffe mixing [5] of p and a\ can be carried out with these parameters, although 
it turns out to be relatively unimportant compared with the BR parametric effect [6] with 
the Dey-Eletsky-Ioffe low-temperature theorem remaining valid. It then becomes obvi- 
ous that the parametric vector (p) mass should be used in the Rapp-Wambach (RW) 
configuration mixing [7] between the "elementary" p (denoted as p e i m ) and the p-sobar 
N*(1520)N~ 1 (denoted as p so bar)- Here the p so f> ar repels the p e i m and they mix, thus mov- 
ing the Psobar down in energy and increasing its p content, while pushing up the p e \ m - In 
other words, BR and RW "fuse" at the lower energies, cooperating in lowering the energy 
of the collective excitation, the p so bar, which goes to zero with chiral restoration [8]. On 
the other hand, BR and RW "defuse", i.e., work in opposite direction for the p eim , where 
BR has to overcome RW in order to move the p-mass down in energy as found empirically 
in the STAR collaboration reconstruction of the p from two emitted pions [6]. Since BR 
scaling and RW configuration mixing give effects of the same general magnitude, it is 
necessary to reconstruct their effects in both the region where they fuse and the one in 
which they defuse in order to disentangle their effects. 

This paper is organized as follows. 

In Section 2 we summarize the highlights of the modern developments. 
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In Section 3 we briefly review the conventional NJL theory, pointing out that the 
cutoff A is sufficiently high to consider constituent quarks but only of the order of the 
p-meson mass, so unsuitable for calculating properties of vector mesons. 

In Section 4 we consider the vector manifestation of chiral symmetry. We show 
that introduction of p-meson degrees of freedom necessitates raising the NJL cutoff by 
a factor of \[2 for consistency in describing the chiral restoration transition. The factor 
arises when p is present because the p-zr loop correction cancels half of the pion loop which 
enters into the running of f n with scale. It follows from the Wilsonian matching that the 
flow of the quark condensate (qq)* governs the movement of the p-mass to zero (in the 
chiral limit). We shall refer to this NJL implemented with the VM as "VM-NJL." 

In Section 5 the Rapp/Wambach configuration mixing is argued, when the para- 
metric structure of the vector manifestation is suitably implemented, to be an additional 
and integral - and not an alternative - ingredient to BR scaling. 

In Section 6 we shall model the modifications in our dropping mass scenario for 
hadronization in RHIC physics. We show that there is a smooth transition with expansion 
of the fireball from gluons into vector and axial-vector mesons, the mesons starting out 
as massless and getting their masses back as temperature and density decrease. We show 
how the interaction of off-shell vector mesons can be reconstructed from lattice gauge 
results for the quark number susceptibility. 

The summary is given in Section 7. 

In the appendices are given some qualitative discussions on recent puzzling obser- 
vations at RHIC in terms of our vector manifestation and gluon-vector-meson "relay" 
scenario. These sections are quite preliminary and will have to be quantified. 

2 The Developments 

The concept of dropping masses in dense/hot medium formulated in 1991 [3] has 
recently undergone rejuvenation in both experimental and theoretical sectors. The original 
idea of implementing the trace anomaly of QCD in terms of a scalar x held which at 
the mean field level gave rise to the scaling of dynamically generated hadron masses 
has recently been reformulated [4] in terms of excitations in dense skyrmion matter, a 
formulation that renders the notion of BR scaling more precise and provides its possible 
link to hidden local symmetry theory with the vector manifestation (HLS/VM) [1]. More 
specifically, the result of [4] provides a connection between the BR scaling and the crucial 
"intrinsic" background dependence brought in by matching to QCD in the HLS/VM 
theory. Furthermore, the x field embedded in medium, i.e., x* \ the average value of 
which represents the radius of the chiral circle and shrinks as chiral symmetry is restored 
is shown by Shuryak and Brown [6] to correctly give the decrease in p-mass with density in 
the pristine atmosphere created in RHIC physics and measured by the STAR collaboration 
just preceding freeze-out of the two pions making up the p-meson. 

This new development is chiefly triggered by a series of work [9] on "vector man- 
ifestation of chiral symmetry" done by Harada and Yamawaki (reviewed in [1]). The 
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work of Harada and Yamawaki is based on the exploitation of hidden local symmetry 
(HLS) that emerges via the gauge-equivalence of nonlinear G/H theory - where G is 
the unbroken global symmetry and H is the subgroup of G that is left unbroken after 
spontaneous symmetry breaking G - to linear G x H\ oca \ theory. This equivalence allows 
them to extend chiral dynamics, i.e., chiral perturbation theory, up through the scale of 
the vector mesons. We might identify this to be a "bottom-up" approach starting with 
a chiral symmetric low-energy theory. How to go up in scale using the gauge-equivalence 
argument is in general not unique; it is known that a variety of gauge theory Lagrangians 
can reduce to a non-gauge theoretic effective Lagrangian, so reversing the reasoning, one 
arrives at the conclusion that a given low-energy effective theory, if unconstrained, can 
flow into a variety of different directions. The crux of the observation made by Harada 
and Yamawaki is that there is a unique flow if the HLS theory is matched a la Wilson 
to QCD at a scale in the vicinity of the chiral scale A x . The matching is made via the 
vector and axial- vector correlators between the hadronic sector described by the effective 
field theory and the QCD sector described by OPE with the gluon and quark condensates 
and color coupling constant given at the matching scale A M . The fixed point structure of 
the effective field theory then locks in the renormalization group (RG) flow towards chiral 
restoration corresponding to the vector manifestation (VM) fixed point 1 . 

We have recently suggested [13, 14, 15] for the three-flavor case that the above 
argument can be made also in a "top-down" theory following Berges-Wetterich's color- 
flavor locking mechanism in the Goldstone mode [16, 17]. If one assumes that quark- 
antiquark pairs condense not only in the color-singlet channel as in the standard picture 
but also in the color-octet channel as was suggested by Berges and Wetterich, then the 
resulting effective theory results in a form identical to the HLS theory of Harada and 
Yamawaki. Since the variables in this theory are the massive quasiquarks and Higgsed 
gluons, it is properly in QCD language and although there is no proof yet, it is plausible 
from the Harada- Yamawaki analysis that the theory possesses the VM fixed point and 
hence the same chiral restoration scenario [15]. 

The vector manifestation viewed both "bottom-up" and "top-down" implies that the 
vector-meson mass goes to zero 2 linearly as (qq) does as the critical point is approached. 
We will refer to this scaling as "Nambu scaling." In fact, in a model calculation connecting 
quantities calculated directly by lattice gauge calculations, Koch and Brown [18] found 
confirmation of Nambu scaling in temperature. These authors considered the increase in 
entropy as 24 hadrons went massless with increasing temperature. (This number of degrees 
of freedom must be preserved in the transition which ends up with a quark multiplicity 

1 It may be helpful to specify how the VM of Harada and Yamawaki differs from Georgi's vector limit 
(VL) [10]. Georgi's VL is characterized by g = 0, a = 1 and /, ^0 whereas the VM is given by g = 0, 
a = 1 and f n = 0. What makes the pion decay constant in the VM flow to zero is the presence of quadratic 
divergence in the renormalization-group equation for / w which will figure importantly in the argument 
we develop below. Unlike Georgi's vector limit, there is no symmetry enhancement in the VM at the 
chiral restoration which is consistent with the symmetries of QCD. It should however be mentioned that 
although perhaps not directly relevant to QCD, Georgi's symmetry enhancement has a significant place in 
the notion of "theory space" being developed for dimensional deconstruction in particle physics [11, 12]. 

2 Unless otherwise noted, we will be mostly considering the chiral limit. 
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of 24 [19]). With (qq)* taken directly from lattice calculations, the entropy S for assumed 
scalings 

M Hadron = ( VjW>_\ ^ Qr (QQ)* ^ ^ 

Mnadron \ (qq) J (qq) 

was compared with the S calculated in the lattice gauge calculation. We show (for com- 
pleteness) the comparison in Fig.l. The rise in lattice gauge entropy through the chiral 
restoration transition is well reproduced with the Nambu scaling. 




Figure 1: Comparison of entropy derived by Koch and Brown from the lattice data on energy 
density (dashed line) and from dropping masses (full lines) calculated in the "Nambu scaling" 
and in the old "Brown-Rho scaling" which is now superseded by a revised one described in this 
paper. In the lower dashed line for (vpip) the bare quark mass has been taken out. Both s and 
(V^V) given as open squares were taken from lattice gauge simulations. 

This linear dependence expected on general ground and observed in the Koch-Brown 
analysis has been verified in the framework of HLS theory only near the critical point both 
in temperature [20] and in density [21]. Recent works [22, 23] indicate that the interplay 
of the "intrinsic temperature/density dependence" 3 and "hadronic thermal/density loop 
effects" makes the situation considerably more subtle than naively thought, so except at 
the two extremes, i.e., at very low temperature/density and very near chiral restoration, 
the scaling behavior is not quantitatively understood in HLS/VM theory. 

3 This effect is missing in all low-energy effective theory calculations so far available in the literature 
except in "unified approach" based on skyrmion soliton model of [4] . 
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It is interesting to note that the strategy involved in the Harada-Yamawaki approach 
is related to that used in QCD sum rule calculations [24]. There the hadronic correlators 
are also matched with the QCD ones at a scale Asr ~ lGeV. A Borel transformation is 
carried out with the result that the scale Asr is replaced by the Borel mass Mb- Thus, 
terms of the form {G ^G^) / M B and {qq) 2 /M B occur in the equations. In parallel to the 
Harada-Yamawaki scheme where the bare parameters determined at the matching scale 
A M are brought down to the on-shell scale m p (or m* in medium) by RGE, the decrease 
in scale in the QCD sum rule approach is then accomplished by bringing the Borel mass 
down to m p (or m*), at which a wide plateau is to be established for the procedure to be 
viable. Indeed such a plateau was found by Adami and Brown in their study of the vector 
meson mass m* in temperature. They found that as {qq)* — > 0, m* — > 0. Thus the zero in 
m* which is a fixed point in HLS/VM, was found in the QCD sum rules by changing only 
the quark condensate, and leaving the color-gauge coupling constant a s and the gluon 
condensate (G^G^) alone. In this formulation, the former changes perturbatively and 
the latter makes "radiative corrections" to the mass generated by the {qq) condensate. A 
crucial question raised in the QCD sum rule calculations of Adami and Brown was why 
Mr should scale with density and/or temperature and whether the assumed scaling in 
order to produce a plateau about m* was correct. We believe this question is answered 
by HLS/VM theory. 

The developments described above clearly indicate that all chiral effective theories 
- regardless of whether they are local gauge symmetric or not - should be so constrained in 
order to correctly represent QCD. In this paper, we argue that the Nambu-Jona-Lasinio 
(NJL) effective theory can be made equivalent to the VM fixed point theory by extending 
it up through the scale of the vector mesons to the scale at which the hadronic theory can 
be matched with QCD and hence can supply the means to interpolate the scaling between 
the two extremes n = and n = n c . To do this, we formulate NJL with constituent quarks 
(or quasiquarks) as variables. We assume that before chiral restoration, transition from 
nucleons to quasiquarks can be made [25, 26] 4 . The phase transition is most naturally 
made in these variables with continuity in the number of degrees of freedom. 

3 Nambu-Jona-Lasinio Theory 

The NJL has the same quadratic divergence in its gap equation, the equation gov- 
erning the dynamical generation of mass, as the VM of Harada and Yamawaki 5 . We wish 
to study chiefly the NJL {qq) condensate, since its medium dependence is found (as we 
discuss below) to run the dropping of the p meson mass. 

Following the work of Hansson et al [27] , we note that in the vacuum, NJL theory 

4 This must be a good assumption as we know the quasiquark picture is fairly good even in free space. 
5 We are using the cutoff to regulate the theory. As discussed by Harada and Yamawaki in [1], the use 
of cutoff does not spoil chiral symmetry in the gap equation while it does in the correlators. 



5 



gives 



N f N c m Q 
An 2 



A4 - vtlq In 



'Al + rriQ 2 
, m Q 2 



(3.1) 



where A4 is the covariant cutoff and tuq is the dynamically generated quark mass of the 
quasiquark. If one takes A 4 = 715 MeV, one obtains, for Nf = 2 and N c = 3, 



{qq) = -(250 MeV) J 



(3.2) 



How the condensate changes in medium in this model has been worked out simply and 
elegantly by Bernard, Meissner and Zahed [28] who found 6 a critical density of ~ 5n 
(where n is the normal nuclear matter density) and T c > 200 MeV with a cutoff of 700 
MeV. We believe this to be ballpark also for the vector manifestation which we discuss in 
the next section. 

In adding nucleons in positive energy states, each added nucleon brings in the same 
scalar condensate as the negative energy nucleon contributes - but negatively - to {qq) 
in eq.(3.1). Thus we get 



(qq)* = ( qq ) + 



6m Q 

47T 2 



kp - m Q 2 In 



kp 



1717 



Qm 



Q 



4tt 2 



A-4 



kp 



m Q 2 ln 



L 4 + m Q z 



A 2 ' — 2> 



kp + m Q < 



(3.3) 



Note that when the Fermi momentum kp reaches the cutoff A4, the condensate vanishes. 
An intuitive interpretation is to view the cutoff as the "Fermi momentum" of the sea of 
the condensed negative energy quarks and that it reflects the magnitude of (qq). The 
condensate vanishes as the positive energy sea becomes as large in magnitude as the 
negative energy sea characterized by A 4 . This would imply a smooth mean field transition 
with the vector mass going to zero with the order parameter. Since the behavior of 
the vector degrees of freedom signals the phase transition, it is likely that there is a 
condensation of these degrees of freedom at the critical density Tl c '^ 1.6., db discontinuity in 
density, the cu-meson being the easiest to condense [29]. We note that (qq) scales as 



(qq) oc m Q k\ 



(3.4) 



for A 4 ^> rriQ, kp, which is the situation up through nuclear matter density no. 

We shall take the usual position on NJL and assume A 4 to be constant, independent 
of density or temperature. This makes it possible to make a connection to the vector 
manifestation of chiral symmetry as we discuss in the next section. 



6 In the chiral limit their critical density was <~ 2no, much too low, but with bare quark mass rh = 
(m u + nrid)/2 w 5 MeV it was more like the 5n that we quote here. See the discussion at the end of 
Section 4. 
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Given A 4 constant, eq.(3.4) shows that in going to finite density, we expect, at low 
density, in the mean field approximation 



m Q oc (qq)* oc /; (3.5) 

with the star representing in-medium quantities. In going from (3.4) to (3.5), we have 
used the fact that tuq in medium satisfies an in-medium gap equation and hence is density 
dependent. The last relation of (3.5) follows from the Goldberger-Treiman relation that 
is expected to hold in medium at mean field in NJL as 7 

G„ QQ m* Q = /* (3.6) 

We have used that the coupling G a QQ does not scale at mean field. We shall refer to the 
scaling (3.5) as "Nambu scaling" [13]. 

In [30] following the argument of Lutz et al [31], it was outlined that from the 
Gell-Mann-Oakes-Renner (GMOR) relation 

fl m l = -m(qq) (3.7) 
where m is the bare quark mass which does not scale with density, 

r:-^- (3. 8 ) 

Although there have been many papers on how the "pion mass" 8 behaves in 
medium, it is a very subtle issue as shown in [32] as the in-medium Goldstone boson 
mass is subject to intricate self-consistency conditions that non-Goldstone bosons are not 
subject to. Given this situation, we take the point of view - supported by [32] - that the 
(parametric) pion mass that appears in current algebra relations is protected (at least) at 
low density by chiral symmetry in which case we would have the scaling 

/x<^n i/2 (3.9) 

that differs from the Nambu scaling (3.5). Empirically, this is much more reasonable at 
low density than the Nambu scaling. This implies that 

rn* p ~ fig* ~ f: ~ (qq)* 1/2 (3.10) 

where we have used that the (hidden) gauge coupling does not scale appreciably for 
n < n . 



7 We are fully aware of the fact that since Lorentz invariance is broken in medium, one has to distinguish 
the time and space components of f v here as well as in what follows. We are keeping Lorentz invariant 
notation but as discussed in [21, 4], the main argument holds valid: This technical matter can be easily 
attended to. 

8 We are being a bit cavalier in what we mean by "pion mass." To be precise, one has to indicate 
whether one is talking about the "parametric mass" or pole mass. Although to the order we consider, 
there is no difference between the two, we have in mind the parametric mass and not the pole mass. 
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On the other hand, Harada, Kim and Rho find that the vector manifestation implies 
close to the chiral restoration density/temperature, 

^~^r Oil) 

m p (qq) 

in agreement with Koch and Brown [18]. What is surprising is that the mass scaling 
seems to accelerate to the form (3.11) at densities only slightly higher than n . As we 
shall see in the next section, m* goes to zero as f*g*yfa, where both /* and g* go to 
zero and a flows towards its fixed point a — 1. In fitting to various theories of nuclear 
matter at densities n > no, Song et al [33, 34] arrived at that the vector (gauge) coupling 
g decreases such that the vector mean field which goes as g* 2 /m v 2 was roughly constant 
above n , although the m v — > evolution is dictated by the flow towards the fixed point. 
This observation combined with that given above suggests that the old BR scaling 

f* m* 

^ « — £ 3.12 



to hold only up to densities where g* begins to decrease appreciably, ^f- dropping from 

then on faster than j-. For the higher densities n > no, where g* scales roughly as to*, 
we find from the flow towards the fixed point that 

<* 9 >* (3.13) 




m p \UJ (qq) ' 



Although this iterative scheme fails if carried to higher order, it does indicate that the 
more rapid Nambu scaling sets in above no- 

From the linear decrease in {qq)* with density [35, 36], 

M! ^ i _ ^ (3.14) 
(qq) flml 

one finds with ct^n = 46 MeV that 

fa)* 



(qq) 

at nuclear matter density, giving 



0.64 (3.15) 



f* 

^-0.8. (3.16) 
We should note that the effective mass of the nucleon m* N /m^r drops faster, by a factor 



9a/ 9a, which is a loop correction going beyond the mean field [3]. (3.16) is consistent 
with the scaling m v /my(no) = $(n ) ~ 0.78 found in the analysis of nuclear matter 
based on the mapping of Landau Fermi-liquid theory to an effective chiral field theory 
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with parameters running a la BR scaling [33, 37, 34, 38]. From our above argument, 
at densities n > n , Nambu scaling sets in, so (/*/ f w ) 2 decreases roughly linearly with 
density from then on, going to zero at ~ 5no- Note that this is roughly the density of 
chiral restoration in NJL quoted above. Thus, the scaling of masses with Nambu scaling 
above no is consistent with the n c obtained in NJL. 

Amusingly, the consequence of this accelerated scaling is that the vector meson mass 
decreases roughly linearly with density n, all the way from n = to n c ~ 5n . 

As n nears n c , eq.(3.3) takes the form 

(qq)* oc m* Q A*A 4 (3.17) 

where A* = A 4 — hp. One can think of A* as an effective cutoff measured relative to the 
Fermi surface that goes to zero as the condensate goes to zero 9 . It is actually a cutoff in 
the phase space available to a positive energy quark before chiral restoration. 

In summary, NJL gives a simple picture of chiral restoration. In mean field, the 
(qq)* condensate goes smoothly to zero, as positive energy quarks are added. Hadron 
masses follow this behavior of the condensate near chiral restoration. This does not 
mean, however, that the phase transition is mean field. There probably is a density 
discontinuity as suggested by Langfeld and Rho due to an induced symmetry breaking of 
vector currents [29]. 



4 The Vector Manifestation of Chiral Symmetry 

Weinberg realized a long time ago that chiral symmetry could be extended to 
straightforwardly include the p-meson in phenomenological Lagrangians [39]. This idea 
was put into a local gauge-invariant form by Bando et al [40] and summarized by Bando, 
Kugo and Yamawaki [41], in their formulation of the nonlinear sigma model as a lin- 
earized hidden local symmetry (HLS) theory in which the p-meson is the gauge particle. 
Since the vector degrees of freedom are strongly coupled to the nucleon, these must be 
included in the hadronic sector up to a scale of A x ~ ^nf-w where they can be matched 
with QCD. The quadratic divergence in Wilsonian renormalization group equations of 
HLS which was recently recognized by Harada and Yamawaki changed the structure of 
the theory radically. In particular it contributes a quadratic term A 2 in the RG flow of 
the order parameter f n and, therefore, is of primary importance in the chiral restoration 
transition. 10 

9 We are grateful to Youngman Kim for pointing out this observation. 

10 Note that power divergences do not appear in chiral perturbation theory in the regime where the 
chiral expansion is valid. However whenever one confronts anomalous scales as when one has quasi- 
bound states or phase transitions, power divergences play the key role. This is the case here in the vector 
manifestation where chiral restoration occurs. Another case that has been discussed in a different context 
is the linear divergence that figures in low-energy nucleon-nucleon scattering in S-wave that has a scale 
invariant fixed point corresponding to an infinite scattering length. See [15] for discussions on this matter 
and references. 
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In a spirit similar to QCD sum rules, a Wilsonian matching is carried out at a scale 
A > 471"/^. This matching determines the bare parameters of the HLS Lagrangian. The 
scale is then systematically brought down all the way to zero by the renormalization group 
equations developed in the effective hadronic sector. The loop corrections are carried out 
with the p-mass counted of the same chiral order as the pion mass. As first pointed out by 
Georgi [10], this is the only way that the p field can be incorporated in a systematic chiral 
perturbation theory. Considering the p-mass as "light" seems unjustified in free space 
but it makes sense when the p-mass drops in medium as mentioned above and developed 
further below. Even in the vacuum, this counting rule combined with the Wilsonian 
matching at A > 47r/ 7r is found to reproduce chiral perturbation theory a la Gasser and 
Leutwyler [42] in irir interactions up to 0(p 4 ). 

The Harada-Yamawaki work is summarized in their Physics Reports [1]. We briefly 
review the parts relevant to our present work. 

Harada and Yamawaki matched the hadronic vector-vector and axial-vector-axial- 
vector correlators defined by 



d 4 xe^(0\TJ^x)J b u (0)\0) = 5 ab {q„q v - g, v q 2 ) U V (Q 2 ), 
d 4 xe^(0\TJ^(x)J b 5u (0)\0) = 5 ab (q,q u - g, u q 2 ) Tl A (Q 2 ) 



with 



Q 2 = -ft 



(4.1) 



(4.2) 



at the scale A > Anf^ between the hadronic sector and the QCD sector. In the hadronic 
sector, they are given by the tree terms 11 



n 



(HLS) 
V 



(Q 2 ) 



F 2 (A) 

Q 2 

F 2 (A) 



M 2 (A) + Q 



-222(A), 

1 - 2g 2 (A)z 3 (A)} - 2zi(A) 



with 



M 2 v (A)=g 2 (A)F 2 (A). 



(4.3) 
(4.4) 

(4.5) 



Here g is the HLS coupling and z^s are 0(p A ) terms in the Lagrangian. In the quark-gluon 
sector, the correlators are given by the OPE to C(l/Q 6 ), 



n 



(QCD) 



n 



(QCD) 
V 



8tt 2 
1 

8^ 



1 + a s /ii) In — + 



Q 2 n 2 { ^ G ^ Gn n3 uogasiqq) 



Q 4 



+ 



3 27 g 6 



i-t , / m Q 2 ,k 2 QG^Gn vr 3 896 a s (qq) 2 
1 + a s /ir)\n— + — - 



Q 4 



3 27 Q 6 



, (4.6) 
(4.7) 



11 Since the correlators are computed at the matching scale, there are no loop corrections. 
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where \x here is the renormalization scale of QCD (e.g., in dimensional regularization). 
Chiral restoration is associated with (qq)* going to zero. In this case Ii^ CD ^ and n|^ CD ' ) 
are equal for all Q 2 . 

The Wilsonian matching equates the correlators of QCD with those of HLS and also 
their derivatives. One obtains 



^( A ) _ rfl d U (QCD) 



A 2 



= -Q' 



dQ 2 



(Q 2 



)2 =A 2. 



(4.8) 



The quark and gluon condensates are obtained phenomenologically as in the QCD sum 
rules, and a s (A) from the perturbative expression. 
The Wilsonian matching implies that 



F 2 M 



a s 27r 2 (^G, u Gn 

7T 3 A 4 



A 2 8tt 2 

which implies in turn that matching with QCD dictates 

f 2 (A) = F n (A) ? 
even though at the point of chiral restoration 

m = o. 



7^0 



(4.9) 



(4.10) 



(4.11) 



Note that i^(A) is a parameter of the Lagrangian defined at the scale A and /^(A) is the 
pion decay constant at the same scale. Thus, /„-(A) = F n (A) is not the order parameter 
but just a parameter of the bare HLS Lagrangian defined at the cutoff A where the 
matching to QCD is made. At chiral restoration (at n = n c or T = T c ) characterized by 
(qq(A; n c ))* = (qq(A; T c ))* = 0, one finds 



5(A) 

o(A) = F 2 (A)/F 2 (A) 
Zl (A)-z 2 (A) 
F.(A) 



0, 

1, 

o, 

0. 



(4.12) 



Here F a is the decay constant for the scalar Goldstone corresponding to the would-be 
longitudinal component of the massive p-meson in the Nambu-Goldstone phase. 

The HLS theory has an elegant fixed-point structure, dictated by the renormaliza- 
tion group equations (RGEs). As in Georgi's vector limit, one has the fixed points 



9 
a 



= 0, 
= 1. 



(4.13) 



Below the p-meson mass the RGE for f n is the same as in the nonlinear chiral Lagrangian 

2N f 



4 fM 



(47T) 



(4.14) 
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where Nf is the number of flavors and [i the renormalization point. The rapid scaling 
with p 2 near the phase transition comes from the quadratic divergences. Chiral restoration 
occurs for 

f«(0) = (4.15) 

at the scale corresponding to pion on-shell. If we integrate (4.14) up to the scale A, say, 
at which the vector meson is integrated out, then 

f = Nf A 2 (4.16) 
U (4vr) 2 { ] 



which gives 



A = ^/.(A). (4.17) 



If we now adopt at this point the KSRF relation m 2 p = 2f 2 gy, then f n (to be evaluated 
at the pion on-shell) is 93 MeV (88 MeV in the chiral limit). Thus 

A = 827(782) MeV. (4.18) 

This is just the cutoff employed in conventional NJL calculations, so everything looks 
fine, except that it is too low a scale to match to QCD and have a description of p-meson 
properties. 

Now once the p-meson is included in the HLS, the (p, 7r)-loop in the axial-vector- 
axial- vector correlator - which is used to calculate f n - cancels half of the pion loop with 
the result that in this theory 12 

A = 4 7 r/ 7r (A). (4.19) 

The NJL with the extended cutoff gives the VM-NJL. Not only is A a factor of \[2 greater 
than A but f n at this higher scale is greater than /„- at m p . Consequently, A > 1 GeV 
is easily employed as the matching scale. Thus, the introduction of the p-meson in the 
loop is necessary to have the cutoff evolve to a sufficiently high scale in order to do the 
Wilsonian matching. 

We should stress that the linear dependence of f n on the cutoff A was already 
contained in NJL, indicating again the relevance of NJL to the Harada-Yamawaki theory. 
The point here is that the cutoff in NJL can be interpreted as effectively playing the role 
of Wilsonian matching scale, as noted before, and the H-Y theory instructs us how to 
extend NJL through m p to 47r/ 7r . 

Since we have the Wilsonian matching, we find that the {(qq)*) 2 terms in U^ CD ^ and 
■q(qcl>) com p} e t e }y dominate the right-hand side, the flow to the fixed point locking their 



12 In order for this cancellation to take place, the scale must be above the p-meson, which is assumed 
in our consideration to be near zero; in other words, the argument here is made in the vicinity of chiral 
restoration. 
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flow toward zero. The other terms in the correlator are down by an order of magnitude. 
Thus, the flow in (qq)* can be obtained from NJL provided the parameters, in particular, 
the cutoff A4, are properly chosen. 

Our strategy is to choose the NJL parameters in order to model results of lattice 
gauge calculations for finite temperature, and then use the NJL results at finite density 
to give the flow to the vector manifestation (VM) fixed point. Although this will be a 
continuing activity to obtain improved NJL parameters as lattice gauge simulations (LGS) 
are improved, we note that it is already possible to close in on such an analysis. The NJL 
calculation of Bernard et al [28] with cutoff A4 = 700 MeV, not far from 715 MeV, gives a 
T c > 200 MeV, essentially that of the LGS of Karsch [43] which is shown in Fig.2. (T c in 
our scheme is unambiguously defined in the chiral limit as the temperature at which the 
hadron masses are zero.) We require that V2A4 > 1 GeV so that the Wilsonian matching 
can be carried out in VM-NJL. 



T (MeV) 

170 210 250 340 510 680 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 

T/T c 



Figure 2: Energy density in units of T 4 for QCD with 3 dynamical light quark flavors [43, 44]. 

We now shortly discuss how the VM-NJL-LGS link can be extended to finite den- 
sities n. For the finite density the explicit chiral symmetry breaking seems important in 
NJL and hence must be taken into account there, because the effects from rh are mag- 
nified by the condensate (qq)* in the Gell-Mann-Oakes-Renner equation (3.8), the (qq)* 
extending well beyond the n c (m = 0). In Bernard et al with fh = 5 MeV, we see that 
(qq)* extends up to n c ~ 5n ; i.e., 

n c (m = 5 MeV) ~ 5n . (4.20) 

We see that even though small bare quark masses enter, they have a major effect on the 
phase transition. This is also true with temperature, where the bare quark masses convert 
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the chiral transition into a smooth crossover. It is thus useful to have VM-NJL since the 
bare quark masses used in the LGS can be put into NJL. 



5 Back to Brown/Rho Plus Rapp/Wambach 

Now that we have established that the p-meson mass must scale with the quark con- 
densate, we need to revisit the Brown/Rho vs. Rapp/Wambach debate in connection with 
the CERES dilepton data. We will propose that the correct approach is "Brown/Rho + 
Rapp/Wambach" with the latter suitably modified to implement the vector manifestation. 

We should recall that the key information from the vector manifestation is that 
when matched a la Wilson to QCD, the bare HLS Lagrangian carries parameters that 
have intrinsic density dependence. Brown-Rho (BR) scaling [3] exploits this density de- 
pendence in formulating it at the constituent quark level [15]. On the other hand, the 
Rapp/Wambach approach [7] is formulated in the language of nucleons and isobars; i.e., 
in the conventional language of nuclear interactions. As such, the Rapp/Wambach is 
immediately accessible through standard nuclear physics with configuration mixing. The 
matching with the experiments is done typically at tree level with phenomenological 
Lagrangians. Now this is of course a well-justified phenomenological approach in very 
low density regimes. But the question arises how this approach can access a higher 
density regime without the account of the intrinsic density dependence in the parame- 
ters of the Lagrangian. In [8], Brown, Li, Rapp, Rho and Wambach suggest how to go 
from Rapp/Wambach to Brown/Rho by endowing the vector meson mass parameter in 
Rapp/Wambach with a suitable density dependence. We reexamine this issue in the light 
of the vector manifestation. 

Rapp et al [45] use a phenomenological Lagrangian of the form 

m p 

to couple the nucleon to the AT* (1520) resonance in the p channel. They use 13 

f P N*N ~ 7 (5.2) 

from fits to the photoproduction data. The N* and N' 1 are combined to have the 
quantum number of the p- meson, the "p-sobar" denoted p so bar- The two branches in the p 
spectral function can be located by solving the real part of the p-meson dispersion relation 
(at q = 0) 

9o = m l + ReX pA r*7v(g )- (5-3) 

Introducing also the backward-going graph, the iV*(1520)A r ~ 1 excitation contributes to 
the self-energy at nuclear matter density n as 

- 3 W^ T { {qo + tTtot /2y-(AEy) (5 ' 4) 
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They actually use f P N*N = 9.2 but their form factor cuts it down to 7 for the relative kinematics. 
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where AE = 1520 - 940 = 580 MeV and T tot = T + T med is the total width of the 
iV*(1520) in medium. When in- medium corrections to the width are ignored, there are 
two solutions to eq.(5.3) 

q~ ~ 540 MeV, g+ ~ 860 MeV. (5.5) 

With inclusion of the widths, both resonances, especially the lower one, are substantially 
broadened. It is the strength in this lower resonance which, because of larger Boltzmann 
factor, gives rise to the excess dileptons at low invariant masses. 

From the Harada and Yamawaki work we learn that the m p in the preceding ex- 
pressions should be replaced by to*, essentially the BR scaling mass 14 . This would 
substantially increase the strength in the lower branch of the p and move it down in 
energy 15 . 

There is considerable disagreement in the literature about the strength of coupling 
of the p to the p so bar- Lutz et al [31] obtain essentially the same value as was obtained 
by Riska and Brown [46] in the chiral quark model f P N*N — 1-44. However Riska 16 
has pointed out that substantial two-quark contributions have been left out in the latter 
work. The work of Nacher et al [47] on the other hand has f P N*N ~ 6, with inclusion of 
form factor for the relative kinematics. Oset suggests in private communication that a 
little stronger p coupling would be even preferred. The work of Nacher et al was used to 
determine how much of the two-pion decay in the 7^ — > 7r + 7r°n process goes through the 
p content of the iV*(1520) excitation in the Langgartner et al [48] experiment. The latter 
authors arrive at a 20% content, in the middle of the 15-25% quoted in the Particle Data 
Physics Booklet. These limits were arrived at from the K-matrix analysis of Longacre 
and Dolbeau [49] and other works quoted in the booklet. Finally, Huber et al [50] find a 
photoproduction cross section for the p° of 10.4 ± 2.5 ph in the three 3 He nucleons in the 
620-700 MeV range of photon energy, relevant for the p so bar, whereas Langgartner et al 
find ~ 3.5 /x6/nucleon for p° production on a proton. In R/W the production per nucleon 
in 3 He should be higher than on a single nucleon; fused with B/R it should be nearly 
double. Helicity measurements enable the Regina group (Huber, Lolos et al) to separate 
out the p° contribution unambiguously and these should remove the possible factor of 2 
uncertainty in the cross section, assuming our theory to be correct. 

It is quite natural that R/W describes most of the "p" participation in lower energy 
nuclear physics. The p so bar comes at zero density at 580 MeV, ~ 3/4 of the elementary 
p mass. It would take a density of ~ |n for BR scaling to bring the elementary p that 
low. With increasing density more and more p is mixed into the p so bar- Quantitatively, 
the amount is increased substantially by the fusing of B/R with R/W. However R/W 
produces too many p-mesons in the region above and about the energy of the elementary 
p, and B/R must be introduced to bring these masses down. 

14 More precisely the intrinsic-density-dependent mass parameter of the HLS Lagrangian M p (A;n) in 
[21]. 

15 Notc that this VM mechanism provides the justification for replacing m p by m* (or more precisely 
the intrinsic mass M*) which was conjectured in [8]. 
16 Private communication from D.O. Riska. 
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6 Hadronization at RHIC 



Based on the Harada-Yamawaki scenario we define T c as the temperature at which 
the vector mesons go massless in the chiral limit. In order to talk about temperature we 
assume equilibration. Hydodynamical descriptions which imply equilibration seem to set 
in fairly early at times r < I fm/c at RHIC [51]. 

As discussed by Brown, Buballa and Rho [25], results of lattice gauge simulations 
indicate that in going upwards through the chiral restoration about half of the trace 
anomaly, amounting to ~ 250 MeV/fm 3 , disappears 17 . This suggests that in going 
downwards through the transition this energy is available to put the hadrons back on 
shell. The feeding back of this energy tends to keep the temperature constant, mimicking 
the mixed phase usually assumed in the phase transition. Recall that the masses change 
with (qq)* near T c , rather than with yj (qq)* relevant at low T. It may take a relatively 
long time for the particles to get their masses back; i.e., to go on shell. In this sense, the 
phase in which the hadrons go back on shell is similar 18 to the long mixed phase [25]. 
This phase may not be, strictly speaking, the mixed phase associated with the first-order 
phase transition but we shall loosely refer to this as such. 

One of the remarkable results from RHIC physics from RHIC physics [52] is the 
common freeze-out temperature of Tf = 174 ± 7 MeV for all hadrons. For a p meson of 
mass m p = 770 MeV this means a total energy of ~ 1090 MeV. 

As developed in [19] quasi-equilibration between n, p and a\ is probably reached 
early in the "mixed phase," before the particles have gone all the way on shell. The p 
continues to equilibrate with two pions p <-> 2n long after the mixed phase, essentially 
down to thermal freeze-out. In fact, because of its large width, Boltzmann factors on 
the two pions distort its shape, essentially lowering its central value by ~ 23 MeV in 
energy and another lowering of ~ 40 MeV comes because it freezes out slightly off-shell 
(due to BR scaling). Because of the above special feature of the p, Braun-Munzinger et 
al [52] do not include it in their Table. However, we believe that the quasi-equilibration 
of the 7r, p and ai forms an underlying framework for the equilibration of the less strongly 
interacting particles. For simplicity let us assume that they all freeze out with a mass 
of ~ 95% of their free space mass in the mixed phase. Then a drop in Tf a of 5 % will 
correct the Braun-Munzinger et all temperature. Given the assumed common scale for 
all masses, we can calculate that they freeze out close to their free-space masses and at 
the Braun-Munzinger et temperature, and then drop all of these by 5 % later. In this 
way, we are doing our calculation as if the particles went back on shell during the mixed 
phase. This applies more particularly to the p-meson. 

At T c each of the massless quark and antiquark coming together to make up the p 
will have the asymptotic energy 3.15T, so the p is formed with energy 6.3T. During the 
mixed phase T does not change, so it must have a mass of 770 MeV at the end of the 

17 These authors attribute this to the Walecka scalar field energy. 

18 We can assess the accuracy of the assumption of constant T in the following way: From the calculation 
of [28] and our argument of Section 3 we see that the particles have about half of their mass back when 
the temperature has dropped 5% from T c . 
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Table 1: Formation and on-shell energies as function of T 



T (MeV) 


£j h in MeV 


6.3T 




(with m p = 770 MeV) 




150 


1070 


945 


175 


1092 


1102 


200 


1145 


1260 



mixed phase where for purposes of calculation we assume it to be on shell. In Table 1 
we list the formation and on-shell energies as function of temperature. We see from the 
table that our assumption of T c ~ Tf Q works well, the two being more or less equal at 
177 MeV. Note that in NJL calculations, T c is essentially unaffected by the bare quark 
mass [28] and our asymptotic value of 6.3T should also be negligibly affected. 

Our view of the phase transition is that the gluons turn into massless vector mesons. 
How this can happen dynamically is outlined in Appendix B which is most certainly not 
the complete story. More theoretically, as proposed in [14], one can think of the gluon- 
meson conversion process as a Higgsing via the color-flavor locking mechanism described 
by Berges and Wetterich [17]. Although the color- flavor- locking transition could very well 
be of first order as argued in [17] for the three-flavor case, in reality with explicitly broken 
chiral symmetry with non-zero quark masses, we see no obstruction to the notion that 
going off and on shell is smooth as we are assuming here. We believe this is consistent 
with the Beg-Shei theorem [53] 19 . In any case, there is a great similarity between massless 
vector mesons and gluons, in that both are confined to the fireball. If one imagines that 
color singlets of gluon pairs go back and forth into pairs of vector mesons, the latter can 
only leave the system when they are back on shell, although they can still scatter on each 
other (with reduced vector coupling g* however) while off-shell in the fireball. We will 
study the interactions of off-shell vector mesons, as they turn into gluons and constituent 
quarks below, using the lattice results for the quark number susceptibility. 

In order to illustrate our principal point made above, namely, that the gluons change 
smoothly into the vector mesons by the Berges- Wetterich type relay, we consider the quark 
number susceptibilities measured in lattice gauge calculations [54], the non-singlet one (in 
the p channel) xns of which is reproduced here in Fig.3, and interpreted by Brown and 
Rho [30]. The description of these in RPA approximation is shown in Fig. 4 in the hadronic 
sector below the chiral restoration point 20 and Fig.5 in the quark/gluon sector above. 
In the analysis of Brown and Rho [30] the vector coupling 21 of the p-meson to nucleons 
at low temperature, g 2 /47r ~ 2.5, gives way to that of constituent quarks as the phase 

19 There are related ideas that come in different context which essentially say the same thing; namely, 
Cheshire Cat Principle, gluon-meson duality etc. This matter is discussed broadly in [15]. 

20 It is interesting to note that this RPA result in NJL - which is a fermionic theory with four-fermi 
interactions - is reproduced by one- loop graphs of the bosonic HLS / VM theory in [23] . 

21 We reserve g for the hadronic vector coupling (e.g., HLS coupling). The color gluon coupling will be 
denoted, whenever needed, as g s , a s etc. 
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Figure 3: The non-singlet (isovector) quark number susceptibility measured by Gottlieb et 
al [54] on lattice. The g appearing here is the color gauge coupling of QCD, denoted in this 
paper as g s . 

transition is approached and then drops to a s « 0.19 (at T = |T C ) after the relay to the 
gluon side. In [30], we conjectured that this reflected the induced flavor symmetry relayed 
to the fundamental color symmetry. This conjecture is supported by the Berges-Wetterich 
scenario [16, 17]. 

Let us study, using Fig.4, how the vector meson goes on shell. We observe that 
the curve for the quark number susceptibility calculated in lattice gauge calculation by 
Gottlieb et al [54] (Fig.3) shows a smooth transition through the critical temperature T c . 
As in Brown and Rho [30] we can parameterize the susceptibility as 

X (T) = (6.1) 

where xo(T) = 2T 2 . For T = 0, we have g 2 /4n = 2.5, so that 

g 2 /m 2 v ^ lOn /m 2 v ^ 31/m^. (6.2) 

Now for T = T c , x ~ 0.5xo — T 2 from which we find 

g* 2 /m v 2 w 2/T 2 w 22/m^ (6.3) 

with my ~ m p and taking T c 220 MeV for the quenched calculation of Gottlieb et al. 
It will turn out that the precise value of T c is unimportant. 
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Figure 4: The quark number susceptibility below the chiral transition temperature is described 
in RPA approximation by summing quasi-quark-quasi-anti-quark bubbles by exchange of p 
mesons or, more properly, by four-Fermi (quark) interactions in the vector-meson channel. 



However, if we continue up to | 



|T C , we see that x ~ -%Xo an d we find 




2 » 
V 



0.2 



(6.4) 




We indicated by the quotation mark an effective constant in the quark-gluon sector that 
can be compared with the quantity in the hadronic sector. We see that there is an order 
of magnitude drop in this quantity in going from the hadronic into the gluonic sector, 
roughly the same as in the g 2 /47r-to-a: s drop. Although the hidden gauge coupling g and 
my go to zero at fixed points at chiral restoration, the ratio is finite and continuous, 
showing that g and my scale in the same way as predicted by the VM. More importantly, 



Figure 5: Perturbative calculation of the quark number susceptibility above the chiral transi- 
tion. The quark and anti-quark coupled to the quark density pq interact by exchanging gluons 
depicted by the wavy line. The gluons are very massive, so the wavy line should be shrunk to 
a point as indicated by the graph with four bubbles. 
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there is a smooth change in u g/my" from the QCD region where the x should be expressed 
in gluons as variables to that in which constituent quarks should be employed as variables. 
The smoothness of this quantity in the lattice gauge simulation of the susceptibility also 
results from the relatively large bare quark masses used in the simulation, giving a smooth 
crossover transition. 

From the dimensional reduction calculation for the gluon side in [30] we see that 
a s ~ 0.19, which grows to g 1 j^n — 2.5 as the gluons are replaced by constituent quarks 
and the p-meson goes back on shell. This may account for the observation that transport 
calculations employing gluons with gluonic couplings are unable to equilibrate sufficiently 
and provide enough pressure for the HBT. 

We remark that we would expect a similar behavior in x as function of density. The 
work of Song et al [33, 34] indicates that from zero density to Uq, the vector gauge coupling 
g remains roughly constant, but my drops to my(n ) / 'my ~ 0.8. However, from then on 
to n c , g*/m v remains more or less unchanged. The decrease in m* without change in g 
was noted some time ago [55] to increase the p-exchange contribution to the tensor force 
in nuclei, thereby reducing the total tensor force, since the p-exchange part has opposite 
sign to that from the pion. This phenomenology supports the scenario that the behavior 
of g with density is similar to that with temperature. 

7 Conclusions 

In this paper, we have shown that the scale of Wilsonian matching of constituent 
quarks to QCD is played by the Nambu-Jona-Lasinio cutoff A 4 . Introduction of the p 
degrees of freedom changes the loop corrections in such a way that the scale is changed 
to y / 2A4 ~ 4:71 f w , the scale used by Harada and Yamawaki in their HLS theory. 

The work of these authors gives a solid basis for BR scaling. Indeed when the VM 
is implemented to NJL, we arrive at an overall picture of how vector-meson masses - 
both parametric and pole - behave from low density to that of chiral restoration that is 
consistent with, but rendered more precise than, the original BR scaling. At very low 
density, chiral perturbation theory is applicable with the p-meson treated as light and at 
densities near chiral restoration, the vector manifestation provides a reliable description. 
In between, our work presents the first indication of how the (BR) scaling behavior inter- 
polates. (This scaling behavior has more recently been reproduced in a description that 
is based on skyrmion structure [4]). At low densities, the parametric F n runs importantly 
while the (HLS) gauge coupling g stays more or less unchanged. So we expect 

m ;/m p ~ fju - my Km)} 1 ' 2 ■ (7.i) 

This means that the phenomenological Lagrangian approach [7] based on Lagrangians 
that have correct symmetries is likely to be qualitatively correct up to ~ Hq. However as 
n exceeds n , the gauge coupling starts falling locked to the VM fixed point and hence 
the phenomenological approach that does not incorporate the flow to the VM fixed point 
cannot correctly access the dynamics in this regime. It is the HLS-type approach matched 
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suitably to QCD that should take over. As one approaches n c , the vector meson (pole) 
mass should fall linearly in {qq)* with the width vanishing as g* 2 ~ ((qq)*) 2 - 22 

We find that for the vector manifestation the density for chiral restoration is n c < 

5n . 

It is found that Rapp/Wambach "widening" should be fused with Brown/Rho scal- 
ing which corresponds to the intrinsic density dependence inherent in the Wilsonian 
matching of the HLS effective theory to QCD. It is likely that this fusing will lower 
n c . In fact, with the coupling between p and [iV*(1520)iv r ~ 1 ] 1 ~ that Rapp and Wambach 
used, n c was estimated to come down to ~ 2.8n [8]. In addition given the difficulty in 
distinguishing between background and subthreshold p contribution, it is clearly too early 
to pin down n c . We can only say at this stage that 2.8 < (n c /n ) < 5 with the upper 
bound being perhaps more relevant because of the explicit chiral symmetry breaking. 

We have reinterpreted our long-standing argument that the transition from the 
hadronic sector to the QCD sector at the chiral transition point must be smooth in terms of 
the Higgsed vector mesons shedding masses as density increases to de-Higgsed (massless) 
gluons and similarly from gapped fermions (baryons) to the gapless quarks etc. There 
is a compelling evidence from the vector isospin susceptibility and axial-vector isospin 
susceptibility that the massless vector mesons (or quasiquarks in NJL) must figure as 
explicit degrees of freedom as do pions near the chiral transition point [23] . The presence 
of such degrees of freedom for temperature near T c makes a completely different prediction 
for the pion decay constants and pion velocity from that of the pion-only scenario. 
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Appendix A: In the Beginning, the Color Glass 

In Section 6 we presented a hadronization scenario based on the Berges-Wetterich 
relay mechanism, in which half of the colored gluons get Higgsed 23 as chiral symmetry 

22 This feature is encoded in the work of Harada, Kim and Rho [21] in the following way. The p-meson 
pole mass is given by the intrinsically running mass M* 2 = a*F* 2 g* 2 plus the dense-loop self-energy term 
g* 2 R where R is a slowly-varying function of density or temperature. Up to n ~ no, g* does not scale, 
so the intrinsic term dominates since the dense-loop corrections are expected to be a small correction. 
Thus the p mass scales like F* which scales as ((qq)*) 1 / 2 - Now for n > no, the intrinsic mass M p scales 
as ((qq)*) 3 / 2 since g* scales as (qq)*. Near the phase transition, therefore, the intrinsic mass term drops 
off fast. Therefore the scaling of the p mass will be dictated by the dense- loop correction which scales 
linearly in (qq)*. 

23 Here we are concerned with the Nf — 2 case of [17]. 
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is broken in the expansion of the fireball formed in the heavy ion collision. 

In a series of papers, McLerran and Venugopalan have developed the color glass 
model, summarized recently in [56]. Here we give a schematic model of how the colored 
glass is shattered, skipping over the detailed picture of McLerran and Venugopalan but 
focusing only on the hadronization scenario. 

With a T ~ 100, RHIC produces matter originally at n ~ 100n . However, only 
one unit of rapidity should be considered, since different units don't communicate with 
each other. Therefore, we have a co-moving n ~ 20no per rapidity interval. Now chiral 
restoration takes place in the time that it takes the nucleon to communicate with the 
condensate which can be described by the exchange of a scalar degree of freedom, say, a 
T xR ~ V(l GeV) = 0.2 fm/c - called U h/Q s " in the color-glass theory, if we include 
the kinetic energy with unit mass of the a. In this time, because of the baryon density, 
the system will be chirally restored, the high baryon density providing the important 
collective effect which distinguishes the nucleus-nucleus collision so markedly from a sum 
of nucleon-nucleon collision. 

The central parts of the colliding sheets of colored glass are shattered through the 
rapid chiral restoration in ~ 0.2 fm/c. Here we will assume the same sort of curve 
for the susceptibility xns as function of density scale as temperature scale: with chiral 
restoration, the couplings will change over from hadronic to perturbative gluonic, thereby 
going to the right side of Fig.2, with small a s gluons. The saturation a s of McLerran- 
Venugopalan is 0.5, somewhat larger than what we arrived at in Section 6. The gluons are 
so copious in the McLerran- Venugopalan theory that each momentum state has many of 
them, and a classical effective field theory can be constructed. Presumably at this stage 
one can talk about temperature because the thermal energy is known to be ~ (E 2 +B 2 )/87r 
assuming equipartition. 

The color glass is shattered into perturbative gluons at about 0.3 fm/c. Quarks are 
less important, being suppressed by order a s as compared with gluons. One can discuss 
the expansion in terms of the entropy current to where r = y/t 2 — z 2 is the proper time, if 
the system is equilibrated. Hydrodynamics [57, 58] shows that in order to obtain sufficient 
flow, equilibration must set in very soon following the color glass phase, at ~ 1 fm/c. The 
subsequent behavior of the hot matter seems to be like that of sticky molasses, rather 
than perturbative quark-gluon plasma. 

The missing part of the RHIC scenario is how equilibration is accomplished from 
the time of the breakup of the color glass down to the hadronization. The large elliptic 
flow t> 2 which develops early tells us that the system does not behave like the predicted 
weakly interacting QGP but more like sticky molasses. This interpretation is supported 
by the low viscosity [59], which implies short mean free path and by the behavior of the 
"balance function" [60] in which the particles born in the decay of a given particle do not 
get very far from each other in rapidity in the time that the fireball evolves. Tantalizingly, 
there are indications from the calculation of hadronic spectral functions about the QCD 
phase transition in the maximum entropy method (MEM) that there are strong collective 
excitations above T c [61]: Strong scalar, pseudoscalar, vector and axial- vector hadronic 
modes are seen. 
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Appendix B: Higgsing of the Gluon; Dynamical Color 
Purifying Filters 

While how the flavor gauge symmetry goes over to the color gauge symmetry at the 
phase transition could be quite complex and intricate depending on the number of fla- 
vors involved, non-zero quark masses and a variety of collective phenomena (e.g., possible 
meson condensation), our scenario is independent of the detailed way in which gluons dis- 
appear with spontaneous break-down of chiral symmetry and the energy goes into mesons 
in the central rapidity (nearly baryon-free) region. The strongly interacting mesons, such 
as the n, p and a x as suggested in [19] will equilibrate as they start to get their masses 
back in going on-shell. 

Although our scenario must happen on general grounds as stressed before, it is 
nonetheless interesting that there are processes that indicate that colored gluons do turn 
into colorless vector mesons. For instance, colliding gluons can be turned into vector 
meson pairs through coupling to Q 2 Q 2 states (where Q stands as before for constituent 
quark) in conventional nuclear physics [62]. As we have argued in the main text that 
the Rapp/Wambach, in conventional nuclear physics language, is only part of the drop- 
ping mass story, the main part being the parametric dependence of the vector meson on 
density/temperature [3, 20, 21], the Li-Liu mechanism that we refer to as "dynamical Hig- 
gsing" may be only part of Higgsing, the main part coming from confinement which forces 
color to disappear with the chiral symmetry breaking (and might be phenomenologically 
mocked up in terms of strongly increasing attractions in color singlet states). 

In the Li-Liu work, two off-shell gluons, one from the pion and the other from the 
proton in a collision such as gives Drell-Yan process, can be converted into two p-mesons, 
via a dynamical Higgsing of gluons, into vector mesons. The two gluons couple through 
Q 2 Q 2 (exotic) states considered a long time ago in [63, 64]. There are a number of these 
states in the energy range > 1 GeV and of particular interest is an / = 2 state at ~ 1.6 
GeV, which is seen experimentally, as well as / = states. The 1 = 2 state can be 
characterized as 

8 

q(Q 2 Q 2 (I = 2)) = aY,V a -V a + bV.V (B.l) 

a=l 

where the V a are colored vector mesons, q r y^\ a q, and V are the usual colorless vector (p) 
mesons. The two gluons couple to * through the first term in the wave function with 
coupling constant a s under the assumption of colored vector dominance adopted by [62] 
and decay into p-mesons through the second term in the wave function. 

This scenario is given credence by the behavior of the 77 — > p°p° cross section [65] , 
to which the gluon-gluon-to-p-p can be directly mapped. Experimentally a sharp peak 
in the cross section is found at E ~ 1.6 GeV, with cross section of a ~ 1CT 7 barns. For 
gluons, a 2 should replace the electromagnetic fine structure constant squared, a 2 , giving 
a factor of ~ 5000 for a s 0.5. Thus, we estimate a ~ 0.5 mb for two gluons to go into 
two p°'s. 
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There must be model states where the V ■ V piece in (B.l) is replaced by n ■ n or 
A ■ A (where A is the axial vector), so there may be substantial direct production of the 
pions and axial vector mesons in the coupled ir, p, a\ system, as well as formation from 
the vector mesons by equilibration. Thus, our "dynamical Higgsing" may have a cross 
section as large as several mbs for gluons of the most favorable energies 

We expect the wave function (B.l) to survive chiral transition without much change. 
The colored fields V a presumably lie at a higher scale and the vector meson mean field 
V ~ g*/my was shown above to increase ~ 25% in going from n = to n = n , and 
then remain constant through chiral restoration. Thus, the Q 2 Q 2 states may, in fact, 
do the dynamical Higgsing of converting gluons into vector mesons, although they must 
be converted ultimately into systems of strongly interacting mesons, with the latter then 
going on shell. 

Incidentally, if the filters are not big enough to purify the color, there may be some 
"sludge" left in terms of glueballs. But these will go out with expansion giving rise to 
decreasing temperature because of small Boltzmann factors. 

Now some caveats. 

Although the group structure in the Higgsing and de-Higgsing can be identified 
unambiguously, it is not feasible to construct the renormalization group flow from per- 
turbative partons to hadrons, although in our scenario the latter would begin massless in 
the chiral limit. The problem is that in practice we are quite far from the chiral limit, 
the constituent quarks having ~ 100 MeV masses in NJL as the dynamically generated 
mass goes to zero. Thus, as seen in the susceptibility, Fig.3, the transition is a smooth 
crossover one. We expect the Higgsing from gluons to vector mesons to take place over 
this crossover. 
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